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Abstract 


The  use  of  superconductors  in  the  field  windings  of 
large  synchronous  machines  displaces  iron  as  a  major 
element  in  the  production  of  the  magnetic  field  and 
promises  smaller,  less  costly  machines  for  a  given  rating. 
The  possibility  of  replacing  the  iron  in  a  transformer 
bank  with  a  rotating  superconducting  field  is  discussed. 
By  placing  the  primary  and  secondary  three  phase  windings 
in  separate  armature  windings  built  in  epoxy  and  fiber- 
glass concentric  shells  around  the  superconducting,  free- 
spinning  field  and  allowing  the  angular  alignment  of  one 
of  tiie  armatures  to  be  varied,  a  unique  transformer  is 
devised.   Unique  features  of  the  synchronous  transformer 
are:   the  ability  to  separately  control  real  and  reactive 
power  flow,  thus  providing  power  factor  correction  while 
simultaneously  acting  as  a  transformer,  and  the  ability  to 
transform  rated  real  power  between  two  power  systems  that 
are  out  of  time  phase,  with  zero  (or  at  a  specified  level 
within  a  range)  reactive  power  flow. 
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INTRODUCTION 

Synchronous  machines  are  commonplace  in  the  modern  world 
and  the  theory  of  their  operation  is  well  developed.   The 
application  of  superconducting  technology  to  the  field  wind- 
ings of  synchronous  machines  is  a  new  and  growing  area  of 
research  [1] .   It  is  expected  that  superconducting  machines 
with  ratings  of  1000  KVA  and  greater  may  have  as  much  as  an 
order  of  magnitude  reduction  in  weight  and  size  over  a  con- 
ventional machine  of  the  same  rating  with  only  about.  0.1 
per  cent  of  machine  rating  required  to  power  the  refrigeration 
unit  of  the  superconductor  [2] . 

Recent  work  at  M.I.T.  anticipates  that  it  may  be  pos- 
sible to  manufacture  and  ship  synchronous  generators  with 
superconducting  field  windings  with  ratings  as  high  as 
10,000  MVA  and  that  construction  and  operation  of  these 
machines  may  be  more  economical  than  current  conventional 
machines  [3] . 

These  economic  advantages  are  brought  about  in  part  by 
the  very  high  flux  densities  with  minimal  power  dissipation 
achievable  with  superconductors,  which  makes  it  possible  to 
eliminate  the  magnetic  circuit  except  possibly  as  an  external 
shield.   Thus  an  expensive  magnetic-iron  rotor  forging  and 
stator  teeth  may  be  eliminated,  which  in  turn  reduces  insula- 
tion requirements  (or  allows  an  increase  in  voltage  rating) 
and  at  the  same  time  increases  the  amount  of  space  that  may 


be  filled  with  armature  conductors.  Elimination  of  most  of 
the  iron  from  the  machine  reduces  per  unit  reactances  which 
may  improve  stability. 

References  [2]  and  [4]  derive  expressions  for  electrical 
parameters  of  synchronous  machines  without  magnetic  shielding. 
References  [5]  and  [6]  extend  these  equations  to  machines 
with  an  exterior  iron  shield. 

In  a  superconducting  alternator  the  superconducting 
field  winding  displaces  iron  as  a  major  element  in  the  pro- 
duction of  the  magnetic  field.   It  might  be  asked  if  this 
can  be  done  for  other  devices,  specifically,  is  it  feasible 
to  replace  the  iron  in  a  transformer  bank  with  a  rotating 
superconducting  field?   Such  a  device,  hereafter  called  a 
synchronous  transformer,  would  have  two  armature  windings 
(one  for  each  of  two  sets  of  three  phase  terminals)  of  the 
form  contemplated  for  a  superconducting  alternator.   It  would 
also  have  a  superconducting  field  winding,  which  is  free 
spinning  as  does  the  field  winding  of  a  synchronous  condenser. 

Electric  energy  is  transferred  from  one  armature  to  the 
other  by  the  interaction  of  the  two  armature  flux  waves  with 
the  field  winding  flux  wave:   one  interaction  being  similar 
to  that  of  a  motor,  the  other  like  a  generator. 

As  it  will  be  shown,  this  proposed  device  has  similari- 
ties with  conventional  transformers,  but  it  has  important 
differences  too.   For  example,  if  the  angular  alignment  of 
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one  of  the  two  armatures  is  made  variable  then  the  device 
has  two  "handles"  on  power  flow:   armature  position  and  main 
field  current.   This  contrasts  with  the  availability  of 
only  one  "handle"  (tap  changing)  with  conventional  trans- 
formers.  This  thesis  will  show  that  these  two  control 
variables  will  enable  the  synchronous  transformer  to  control 
separately  both  real  and  reactive  power  flow.   Thus  it  can 
act  as  a  power  factor  correction  device  while  simultaneously 
acting  as  a  transformer.   Additionally,  it  can  transform 
rated  real  power  between  two  power  system  buses  that  are  out 
of  time  phase  with  no  reactive  power  flow  or  within  a  range 
of  desired  reactive  power  levels. 

Such  a  device  could  be  constructed  by  embedding  the 
armature  windings  in  epoxy  and  fiberglass  to  a  fiberglass 
(or  other  non-magnetic  material)  pipe  for  structural  strength 
[7,6].   Two  such  armature  windings  would  be  required.   The 
axis  of  the  outer  armature  would  be  adjustable  so  that  it 
could  be  rotated  about  the  longitudinal  axis  of  the  rotor 
(and  thus  also  about  the  longitudinal  axis  of  the  inner 
armature)  by  an  external  control,  but  it  would  not  be  free 
to  spin.   The  returning  flux  path  and  shielding  would  be 
provided  by  a  concentric  iron  shell  external  to  the  two 
armature  shells.   Figure  1  illustrates  the  physical  configu- 
ration of  the  synchronous  transformer.   Figure  2  illustrates 
the  electrical  configuration. 
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Figure  1 

Physical  Configuration  of  the  Synchronous  Transformer/ 
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Figure  2 

Electrical  Configuration 
Armature  2  is  Shown  Rotated  by  an  Angle  <J>  from  Armature  1 


12 

There  is  a  growing  interest  in  the  use  of  supercon- 
ducting machinery  for  marine  propulsion.   In  19  66  Sibley 
et  al  [8]  estimated  some  potential  advantages  of  super- 
conducting electrical  machinery  over  conventional  propulsion 
plants.   In  a  later  work  Laine  [9]  considered  the  possibil- 
ity of  using  a  superconducting  transmission  link  in  a 
destroyer-type  ship.   Matjasko  [10]  investigated  a  control 
system  for  a  destroyer  propelled  by  superconducting  electric 
machinery.   Greene  [11]  examined  transmission  systems 
including  a  cycloconverter  scheme  between  a  synchronous 
generator  and  an  18,000  SHP  (13.5  MW)  synchronous  motor 
driving  a  propeller. 

While  in  its  present  configuration  it  would  seem  that 
the  synchronous  transformer  would  perhaps  be  more  useful  to 
the  electric  power  industry  than  to  marine  propulsion,  it  is, 
however,  similar  in  configuration  to  a  class  of  dual  arma- 
ture superconducting  machines  which  may  have  great  potential 
as  an  electric  propulsion  ship  system.   This  class  of  AC 
machines  may  have  the  control  characteristics  of  a  DC  machine 
system  but  with  the  additional  advantage  that  a  solid  state 
switching  circuit  will  replace  the  bulky  cycloconverter. 
This  class  of  dual  armature  machines  will  not  be  discussed 
further  here  as  they  are  the  proprietary  information  of  the 
Cryogenic  Engineering  Laboratory  at  M.I.T. 


Table  1  13 

GLOSSARY  OF  TERMS 

Note:  Unless  otherwise  noted,  the  following  subscripts 
will  serve  as  modifiers  to  the  symbols. 


Subscripts 

a  phase  a,  armature  1 

b  phase  b,  armature  1 

c  phase  c ,  armature  1 

a  phase  a,  armature  2 

8  phase  6 ,  armature  2 

y  phase  y ,    armature  2 

d  direct  axis,  armature  1 

q  quadrature  axis,  armature  1 

6  direct  axis,  armature  2 

£  quadrature  axis,  armature  2 

f  field 

o  zero  sequence,  armature  1 

ft  zero  sequence,  armature  2 

P,  phase  variables,  armature  1 

P-  pnase  variables,  armature  2 

R,  rotor  variables ,  armature  1 

R_  rotor  variables,  armature  2 

L,  armature  1  to  field 

L_  armature  2  to  field 

t,  equ:  valent  at  the  terminals  ,  armature  1 

t0  equivalent  at  the  terminals ,  armature  2 


Glossary  of  Terms 
(continued) 
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Symbo 

Is 

i 

current 

V 

voltage 

A 

winding  flux 

L 

inductance 

e 

angular  posi 

angular  position  of  the  rotor  from  the  axis  of 
phase  a,  armature  1,  positive  in  direction  of 
rotation 

cf)      angular  position  of  the  axis  of  phase  a, 

armature  2,  from  the  axis  of  phase  a,  arma- 
ture 1,  positive  in  the  direction  of  rotation 

M     maximum  value  of  mutual  inductance  between 

phase  a,  armature  1,  and  phase  a,  armature  2 

(i.e.,  L    at  d>  =  0) 
aa     r 

L,      maximum  value  of  mutual  inductance  between 

phase  a,  armature  1,  and  the  field 

(i.e.,  L  £  at  8  =  0 
af 

L~     maximum  value  of  mutual  inductance  between 

phase  a,  armature  2,  and  the  field 

(i.e.,  L  r.    at  0  =  <£) 
a  l. 

6       spatial  angle  from  the  rotor  axis  to  the  arma- 
ture flux  wave  of  armature  1,  positive  in  the 
direction  of  rotation. 

6       spatial  angle  from  the  armature  flux  wave  of 
armature  1  to  the  armature  flux  wave  of 
armature  2,  positive  in  the  direction  of 
rotation 

V,  terminal  voltage  of  armature  1  in  per  unit 

V«  terminal  voltage  of  armature  2  in  per  unit 

x  per  unit  reactance 

E  air  gap  voltage 
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RESULTS  OF  A  STATIC  ANALYSIS 

For  any  winding  of  a  synchronous  machine  the  voltage  at 
the  terminals  is  equal  to  the  product  of  the  winding  resis- 
tance and  current  through  the  winding  plus  the  time  rate  of 
change  of  the  winding  flux.   For  any  winding  the  flux  is  the 
summation  of  the  self   inductance  times  the  current  through 
the  winding  and  the  product  of  mutual   inductances  with 
other  windings  (which  may  depend  on  rotor  angular  position) 
and  the  currents  through  those  windings. 

Thus: 

dA 
Vm  =  r  i   +  -rr^  (1) 

m    mm   dt 

A   =  7  X   i  (2) 

m    L      nm  n 
n 

Losses  in  large  synchronous  machines  are  less  than  one 
per  cent  of  machine  rating  and  the  assumption  of  a  lossless 
model  does  not  have  much  effect  on  the  static  analysis. 
Therefore  one  can  assume  that  winding  resistances  are 
neglible  yielding: 

dA 
Vm  "  diT  (3> 

The  machine  flux  current  relationships  may  be  expressed 
by: 
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a    ab    ac    aa    a3    ay    af 


Lab   Lb    Lbc   Lba   ^3   LbY 


bf 


ac   be    c     ca    c3    cy    cf 


aa   ba    ca    a     a3    ay    af 


LaB   Lb3  Lc3   La3   L3    L3y   L3f 


L     L,     L     L     L„    L     L  - 
ay   by   cy   ay    3y   y    yf 


Laf  Lbf  Lcf   Laf   L3f  Lyf   Lf 


a 


H 


y 


(4) 


Because  of  fixed  geometry  the  following  relations  apply: 


L   =  L,  =  L 
a    b     c 


(5) 


L   =  L0  =  L 

a    3    y 


L  ,  =  L    =  L, 

ab     ac     be 


(6) 
(7) 


a3     3y     ay 


(8) 


It  is  important  to  note  that  the  geometry  between  arma- 
tures is  not  fixed  since  armature  2  may  be  rotated  and  since 
phase  axes  are  120°  apart. 


aa     a3     ay 


(9) 


It  is  assumed  that  all  rotor  to  armature  mutual 
inductances  vary  with  the  cosine  of  rotor  position, with  the 

maximum  values  of  L  ,  and  L  -  when  the  rotor  axis  aligns  with 

af       af  r 

the  respective  phase  axes  of  armature  1  and  armature  2.   This 
is  the  same  as  neglecting  higher  order  space  harmonics  in  the 
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rotor  air  gap  flux  and  assuming  a  non-salient  rotor. 
It  is  further  assumed  that  armature  2  to  armature  1  mutual 
inductances  vary  as  the  cosine  of  the  relative  angular  posi- 
tion of  the  respective  phase  axes  with  the  maximum  value  of 
L   when  the  armature  displacement  angle,  4> ,  equals  zero. 
This  value  L    (maximum)  will  be  denoted  by  M.   This  assump- 
tion is  equivalent  to  neglecting  higher  order  space  harmonics 
in  the  air  gap  flux  between  armatures.   Higher  order  space 
harmonics  die  away  more  rapidly  than  the  fundamental.   Since 
this  is  an  air  gap  machine,  disregarding  higher  order  space 
harmonics  in  the  air  gaps  will  introduce  only  a  small  error 
[6]  .   Substitution  into  equation  (4)  yields  equation  (10) . 

Equation  (3)  gives: 


V 


V, 


V 


V 


V, 


V 


V, 


_d_ 
dt 


a 


(11) 


Y 


At  this  point  it  is  convenient  to  transform  equation 
(11)  into  the  rotating  coordinate  system  of  the  rotor.   The 
subscripts  d  and  q  correspond  to  quantities  of  stator  1  , 
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rotating  with  the  direct  and  quadrature  axis  of  the  field. 
The  subscript,  o,  corresponds  to  zero  sequence  quantities  of 
stator  1,  that  is  in  phase  in  all  three  windings  of  stator  1. 
The  corresponding  subscripts  for  stator  2  are  6,  E, ,    and  ft. 
The  details  of  this  transformation  are  given  in  Appendix  A. 

The  result  of  transforming  to  rotor  coordinates  yields 
equations  (12)  and  (13). 

By  comparison  of  equations  (12)  and  (13)  the  following 
relationships  can  be  recognized: 

Vd  =  £  (Xd>  -  UXq  (14) 

Vq  =  E  (Y  +  "Xd  <15) 

V6  =  St     (V  "  uX?  (16) 

v?  =  £    (X5)  +  UX6  (17) 

vf  -  re  xf  (18) 

It  can  be  seen  that  zero  sequence  variables  can  be 
eliminated  from  equations  (12)  and  (13)  since  they  are  com- 
pletely uncoupled  from  other  variables. 

For  ease  in  notation  the  following  definitions  are 
introduced: 

Lw  =  L   -  L  .  (19) 

E     a     ab 

L   =  L   -  L  r  (20) 

£    a     a£ 

Thus  equation  (13)  becomes: 


0 


0 
L, 


3  3 

-=■  M   cos <J)      -  ■=■  M   sin4> 


3  .  3 

■j  M   sine})      ■=■  M   cosif) 


3  3 

y  M   cost})      2"  M   sintj)      L 


3  .  3 

•j  M    sincj)      -J  M   costj)       0 


0 

L 


/  2  Ll 


71  L 
/  2  L2 


/  2  Ll 


71  L 
/  2  L2 


0 
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x5 


(21) 


Now  that  the  basic  relations  governing  the  synchronous 
transformer  have  been  developed,  the  operation  of  the  trans- 
former between  two  buses  will  be  looked  at.   Motor  sign  con- 
vention will  be  observed  throughout,  that  is,  real  power 
input  to  the  machine  from  either  system  (bus)  is  defined  as 
positive  and  lagging  current  determines  the  positive  direc- 
tion for  reactive  power. 

Since  the  flux-current  relationships  have  been  trans- 
formed into  rotor  d-q  coordinates,  the  following  definitions 
of  armature  equivalent  flux  vectors  are  introduced: 


Xtl  =  /Xd'  +  V 


Xt2  "  ^V  +  V 


(22) 

(23) 


Similarly,  the  equivalent  armature  voltage  vectors  are 
likewise  defined: 


V_  =  /v,z  +  V  2 
tl      d      q 


Vt2  =  /V6'  +  V 


(24) 
(25) 
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The  angle  5   is  defined  as  the  spatial  angle  from  the 
a 

rotor  axis  to  the  armature  flux  of  armature  1,  positive  in 

the  direction  of  rotation.   The  angle  <5   is  defined  as  the 

o 

spatial  angle  from  the  armature  flux  wave  of  armature  1  to 
the  armature  flux  wave  of  armature  2,  positive  in  the  direc- 
tion of  rotation.   Figure  3  illustrates  these  sign  conven- 
tions . 

Alternatively,  since  the  voltage  vectors  are  a  ninety 
degree  rotation  of  the  armature  flux  wave  vectors,  these  same 
definitions  can  be  viewed  as  applying  equally  between  the  air 
gap  voltage,  terminal  voltage  of  armature  1  and  terminal 
voltage  of  armature  2.   Figure  4  illustrates  the  sign  conven- 
tions of  angles  6   and  5   in  terms  of  the  voltage  vectors. 

The  following  definitions  follow  directly: 


V ,  =  -  V, ,  sin6=  (26) 

d      tl      a 


V  =  V\,  cos6^  (27) 

q     tl      a 


V.  =-V.0  sin(5   +  6  )  (28) 

o     t2       a     o 


Vr  =  V.~  cos(6=  +  6  )  •   ,(29) 

£,     tz       a     o 


If  the  restriction  to  strict  steady  state  operation  is 
assumed,  then  time  derivative  terms  become  zero.  Equations 
(14)  through  (17)  become: 
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Figure  3 

Illustration  of  Sign  Conventions 

for  Spatial  Angles  5   and  6 
^        ^      a      o 
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Figure  4 

Illustration  of  Sign  Conventions  for  Spatial  Angles 

6   and  5  in  Terms  of  Voltage  Vectors 
a     o  ^ 
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Vd   =    -    u)Aq  (30) 


Vq   =    WXd  ■  (31) 


Vg    =    -    coX  (32) 


V      =    o>A6  (33) 


Substituting   of    equations    (21)    through    (29)    into 
equations    (30)    through    (33)    yields: 

3  3 

Vtl    cos<5a   =    ~c0^LEi      +    7   M   sin<i):i-(5    +    j   M   cos^ir)  (34) 

3  3  /3~ 

Vfcl    sin(5a   =    UJ(LEid    +    2   M   cos<f>ig    ~   T  M   sin^ir    +    /  J  Llif^     ^35^ 

V,  0    cos(6,    +    6    )    =    -a)(-   4  M   sin<|>i,    +  1m   coscbi      +   L    ir)        (36) 
t.z  ao  z  a        z  q  £    s 

3  3  /3~ 

V.  -    sin(6    +6    )    =    to(^  M   cos(J>i,    +   ^  M   sin4>i      +   L.  i.    +     /  ■=■  L„i,. 
t.Z  ao  2  a         2  q  £    o  /    2       2    f 

(37) 


Since  the  rotor  is  free  spinning  there  will  be  no 
electrical-mechanical  power  transfer  from  the  rotor.   If  we 
ignore  all  losses  (friction,  windage,  copper  losses,  etc.) 
then  the  net  power  of  the  synchronous  transformer  will  be 
zero. 


0  =  V,i,  +  V  i   +  v.!.  +  Vrir  (38) 

dd     qq     o    6  t,    E, 


This  is  not  to  say  that  the  device  cannot  absorb  power 
from  one  bus  under  the  lossless  assumption,  but  rather 


equation  (38)  tells  us  that  real  power  into  either  _ 

armature  must  go  out  the  other  armature. 

Substitution  of  equations ( (14)  through  (17)  with  the 
steady  state  asumption  d/dt  =  0  and  equation  (21)  into 
equation  (38)  yields: 


0  =  wLji   +  wL2i  (39) 


Equations  (34)  through  (37)  and  equation  (39)  are  the 

governing  equations  for  the  synchronous  transformer  operating 

between  two  buses.   Independent  variables  are  if,  4> ,  V, 

and  V„.   Dependent  variables  are  i,,  ixi    i  ,  ir,  6    ,  and  <$  . 
2      ^  d    &    q    €   a   .    o 

Thus  we  have  five  equations  and  six  unknowns.   The  method  of 

solution  used  is  to  assign  values  of  6   for  a  sufficiently 

o  u 

wide  range  and  at  each  value  of  5   solve  the  resulting  equa- 
tions whicn  then  have  only  five  unknowns.   However,  before 
proceeding  to  the  solutions  of  these  equations  it  is  con- 
venient to  put  them  into  the  per  unit  system.   Details  of 
this  are  in  Appendix  B. 

The  governing  equations  for  the  synchronous  transformer 
operating  between  two  buses  in  per  unit  are: 


V.  sin6  =  x,i   +  x  i,  sin6  +  x  i_  cosd>  (40) 

1      a     dq     m6  m  ^ 

V,  cos6  =  x.,i,  +  x  ir  coscb  -  x  i_  sind)  +  xxli_    (41) 
1      a     dd     m6     r     m^  Llf 

V^  sin(6  +6  )  =  -  x  i,  sind)  +  x  i   cosd)  +  xrir  (42) 

2       ao        md     r     mq  6  E, 

V»  cos  (6  +5  )  =  x  i  ,  cosd)  +  x  i   sind)  +  xa ,  +  xT  _i£    (43) 
2       ao      md     r     mq  6    6  L2  f 

0  =  xT  i   +  x   i  (44) 

LI  q     L2  t, 
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In  order  to  use  Gauss-Seidel  computer  aided  iterative 

solution  these  equations  were  rewritten  as  follows: 

x.,i   +  x  ir  sincb  +  x  i^  cosd) 

6   =  sin  x  l  V,  J  ^q*} 

a  1 


Vl  CQs6a  -  V-6  COS^  +  Vg  sin*  -  xLlif 
id — (46) 


V~  sin(6  +6  )  +  x  i  ,  sind>  -  x  i   cosd) 
2       a   o      md     r     mq     r  .  .  _, 

H  = ^ (47> 


V2  cos(ia+6o)-  xmid  cosj,  -  xmiq  sinf  -  xL2if    (48) 
x6  x6 


q      XL1 


Appendix  (C)  discussed  the  Gauss-Seidel  iterative  solu- 
tion.  Reference  12  discusses  this  and  other  applicable 
methods  in  great  detail. 

Reasonable  values  of  per-unit  inductances  for  a  syn- 
chronous machine  with  a  superconducting  field  winding  were 
chosen  [6,  13].   The  assigned  inductances  are  listed  in 
Table  2. 
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Table  2 

PER-UNIT  REACTANCES  USED  IN  STATIC  ANALYSIS 


Armature  1  Synchronous  Reactance      x,  =  0.5  p.u. 

Armature  2  Synchronous  Reactance      x-  =  0.5  p.u. 

Mutual  Reactance  Between  Armatures     x  =  0.4  p.u. 

m  c 

Armature  1  to  Field  Resistance       x  ,  =  1.0  p.u. 

Li  J. 

Armature  2  to  Field  Reactance         x  0  =  1.0  p.u. 


A  computer  program  was  prepared  to  solve  simultaneous 
equations  (45)  through  (49)  by  the  Gauss-Seidel  iterative 
method  for  a  range  of  armature  displacement  angles,  <£  ,  from 
-9  0°  to  +9  0°  and  a  range  of  spatial  angles  between  armature 
voltage  vectors,  6  ,  from  0°  to  90°  and  for  different  levels 
of  field  current.   In  all  cases  reactances  from  Table  2  and 
one  per  unit  terminal  voltages  were  used.   For  each  level  of 
field  current  and  each  choice  of  angles  <J)  and  6   the  computer 
calculated  and  printed  values  of  direct  and  quadratic  axis 
currents  and  voltages  for  both  armatures,  real  and  reactive 
power  for  both  armatures ,  and  the  spatial  angle  between  the 
air  gap  voltage  vector  and  the  voltage  vector  of  armature  1. 
This  produced  a  wealth  of  data  from  which  Figures  5  through 
10  were  extracted. 

It  is  important  to  note  here  that  while  the  angles  &a 
and  6   are  spatial  angles  relating  the  voltage  vectors  to  the 
rotor,  the  ang'lf  (J)  is  a  physical  angle  by  which  the  armatures 
are  displaced..  Thus  the  time  phase  angle  between  the  power 
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systems  attached  to  the  armatures  is  a  combination  of  their 
physical  and  electrical  angular  difference,  6   +  <p.      Since 
we  are  not  interested  in  internal  angles  to  characterize  the 
device  this  power  system  time  phase  angle  will  be  the  sig- 
nificant angle  used  to  characterize  the  device. 

Since  we  are  primarily  interested  in  the  synchronous 
transformer  to  transfer  power,  the  solutions  yielding  real 
power  equal  to  one  per  unit  were  collated  and  graphed  in 
Figures  5  through  7.   Figures  8  through  10  were  constructed 
from  Figures  6  and  7.   Figure  8  was  obtained  by  cross  plot- 
ting the  armature  displacement  angle,  cf> ,  from  Figure  6 
corresponding  to  the  zero  crossings  of  reactive  power  flow 
from  Figure  7.   Thus  Figure  8  gives  the  setting  points  for 
armature  displacement  angle,  c}>,  and  field  current,  If,  for  a 
range  of  power  system  time  phase  differences  in  which  one  per 
unit  real  power  can  be  transferred  with  no  reactive  power. 
Figures  9  and  10  yield  similar  results  for  reactive  power 
flow  equal  to  plus  and  minus  0.1  per  unit. 
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Figure  5 

Curves  of  armature  displacement  angle,  cj> ,  versus 
spatial  angle  between  armature  voltage  vectors  for 
one  per  unit  real  power  transfer  and  various  levels 
of  field  current,  1~ 
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Figure  6 

Curves  of  armature  displacement  angle,  <j> ,  versus 
power  system  time  phase  angle  difference,  60  +  <J> , 
for  one  per  unit  real  power  transfer  and  various 
levels  of  field  current,  1^ 
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Figure  7 


Curves  of  reactive  power,  Q,  versus  power  system 
time  phase  angle  difference,  6   +  $ ,  for  one  per 
unit  real  power  transfer  and  various  levels  of 
field  current,  I- 
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Figure  8 


Armature  displacement  angle,  <J> ,  and  field  current, 
If,  versus  power  system  time  phase  angle  differ- 
ence, 6Q    +  4>,  for  one  per  unit  real  power  transfer 
and  reactive  power  equal  to  zero. 
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Figure    9 

Armature  displacement  angle,  <b ,  and  field  current, 
If,  versus  power  system  time  phase  angle  differ- 
ence, 60  +  $,  for  one  per  unit  real  power  transfer 
and  reactive  power  equal  to  +  0.1  per  unit. 
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Figure  10 

Armature  displacement  angle,  <£ ,  and  field  current, 
I-,  versus  power  system  time  phase  angle  differ- 
ence, 60  +  <j> ,  for  one  per  unit  real  power  transfer 
and  reactive  power  equal  to  -0.1  per  unit. 
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DISCUSSION  OF  RESULTS 

The  important  result  of  this  thesis  is  that  the 
synchronous  transformer  can  be  characterized  by  a  set  of 
curves  in  which  the  armature  displacement  angle  and  field 
current  can  be  determined  for  a  given  power  system  time 
phase  angle  difference  for  operation  at  real  power  equal  to 
one  per  unit  and  at  various  levels  of  reactive  power.   In  an 
ordinary  transformer  there  is  always  reactive  power  flow  due 
to  leakage  reactance,  but  in  the  synchronous  transformer 
settings  may  readily  be  chosen  for  unity  power  factor  power 
transfer.   Additionally,  the  synchronous  transformer  can 
operate  at  unity  power  factor  when  transferring  power  between 
two  systems  that  are  out  of  time  phase  or  it  can  simul- 
taneously act  as  a  transformer  and  a  synchronous  condenser 
between  two  buses  that  may  be  in  or  out  of  time  phase.   Thus 
it  is  seen  that  a  synchronous  transformer  in  addition  to 
being  able  to  act  as  a  transformer  and  as  a  synchronous 
condenser,  can  also  act  as  a  phase  shifter. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

The  results  of  this  static  analysis  have  been  encourag- 
ing.  It  has  been  shown  that  the  analytic  model  of  equations 
(40)  through  (44)  yield  a  set  of  curves  in  which  the  armature 
displacement  angle  and  field  current  can  be  determined  for  a 
given  power  system  time  phase  angle  difference  for  operation 
at  real  power  equal  to  one  per  unit  and  at  various  levels 
of  reactive  power.   The  synchronous  transformer  can  transfer 
real  power  from  one  bus  to  another  whether  its  time  phase 
is  behind,  equal  to,  or  ahead  of  that  of  the  other  bus. 

Recommendations 

It  is  recommended  that  further  investigations  of  the 
synchronous  transformer  be  pursued.   Further  efforts  should 
include  a  dynamic  analysis  with  faults  and  transients  and 
an  economic  analysis  to  determine  the  best  level  of  its 
usage  within  the  power  systems  and  its  general  economic 
feasibility. 
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APPENDIX   A 
TRANSFORMATION    TO    ROTOR   COORDINATES 


TRANSFORMATION  TO  ROTOR  COORDINATES 
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The  machine  voltage  current  relationships  (equations 
(10)  and  (11)  may  be  expressed  by  equation  (A-l)  (page  43) . 

Define  a  pair  of  reciprocal  transforms  relating  phase 
and  rotor  corrdinates ,  where  u  represents  a  current  or  a 
voltage.   These  transforms  are  similar  to  those  of  Lewis  [14] 


u 


u 


u 


o  -J 


=       T 


u 


u, 


u 
L      c 


(A-2) 


u 


u. 


u 


=    T 


-1 


u 


u 


u 


(A-3) 


T  =    /£ 

T    "    /    3 


cosG       .cos(6    -    2tt/3)  cos(0   +    2tt/3) 

-sine      -sin(6    -    2tt/3)       -sin(6    +    2tt/3) 
1/2  1/2  1/2 


(A-4) 


T    E 


/2 
/   3 


cose  -sine  1 

cos(6    -    2tt/3)       -sin(6    -    2tt/3)       1 
cos(6    +    2tt/3)       -sin(6   +    2tt/3)       1 


(A-5) 


By  partitioning  equation  (A-l)  as  shown  we  can  reduce 
the  notation  as  shown  below  and  treat  the  partitions 
separately. 
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PI 
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P2 


V, 


_d_ 

dt 


A,    A     A  - 
a     aa    af 


A     A     A  c 
aa    a     af 


A.    A,    A. 

fa    fa    f 


'PI 


P2 
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(A-6) 


Note  that  these  submatrices  are  not  symmetric,  i.e., 
A    ^  A^  .   The  off  diagonal  elements  are  transposes  of  their 

a  a. 

counterparts. 

The  transformed  version  of  equation  (A-6)  is: 


V 


Rl 


V 


R2 


V, 


T-^A  T~l         T-j^A   T  1       T-£-  A  _ 
dt  a        dt  aa       dt   af 


T-^A   T  X   T~  AT1   T-£-  A  _ 
dt  aa       dt   a       dt   af 


dt  AfaT     dt  AfaT     dt  Af 


■Rl 


R2 


(A-7) 


Carrying  out  the  above  matrix  operations  yields  equa- 
tion (12)  of  the  text. 

If  the  flux-current  relationships  (equation  (10))  are 
partitioned  in  the  same  manner  as  equation  (A-l) ,  they 
become : 


PI 


P2 


A      A  - 
aa     af 


A      A 

aa    a 


A 


af 


fa     fa     f 


PI 


"P2 


(A- 8) 


The  transformed  version  of  equation  (A-8)  is 
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•Rl 


R2 


TAT 
a 


-1 


T   A         T 

aa 


-1 


T    A         T 

aa 


T   A         T 

a 


-1 


-1 


T    A 


af 
T   A 


af 


A    _       T 
fa 


-1 


A£      T 
fa 


-1 


Rl 


■R2 


(A-9) 


Carrying   out    this    matrix  operation  yields    equation 
(13)    of   the  text. 
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APPENDIX  B 
THE  PER-UNIT  EQUATIONS 


£7_ 

THE  PER-UNIT  EQUATIONS 


There  are  various  techniques  for  changing  equations  into 
per-unit  equations,  all  of  them  related.   The  method  that  I 
shall  employ  parallels  that  of  the  notes  for  Power  Systems 
Engineering  II,  course  6.552  at  M.I.T.,  by  Prof.  G.  L. 
Wilson. 

The  following  base  quantities  are  specified: 

Axis    Base  Voltage    Base  Current 


d 

Vo 

I 
o 

q 

Vo 

N 

I 
o 

N 

6 

V 

a 

I 

a 

o 

N 
a 

N 

o 

N 
a 

N 

£ 

V 

a 

I 

a 

o 

N 

o 

N 

a 

a 

N^ 

N 

f 

3 

a 

f 

V 

I 

o 

N~ 
a 

2 

o 

*f 

The  subscripted  N's  in  the  above  base  quantities 
represent  turns  ratios. 

The  equations  to  be  put  into  per-unit  are  the  governing 
equations  of  the  synchronous  transformer  operating  between 
two  buses,  equations  (34)  through  (37)  and  (39).   These  are 
listed  below: 
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3  3 

V1    cos    6a  =    -co  (Lgi      +   ^   sin*    i6    +    9^   cos4>    i-r)  (B-l) 

V1   sin   6a  =   w(LEid  +   |m  cosd)   i6    -   |m  sin<j>    i£  +    /f  L1if) 

(B-2) 

V0    cos  (5    +6    )    =  -w(-   |m  sincj)    i,   +   |m   cost})    i      +   L    i_)  (B-3) 

2  ao  z  d         2  q  e    £ 

V0    sin(5    +6    )    =    oj(|m   cost})    i,    +   |m   sine})    i      +   L    i-    +     f\  L-iJ 
2  ao  2  d2  q  e    6         /    2      2    f 

(B-4) 

0      =      coL,i      +    coL^i  (B-5) 

1    q  2Z 


However,  we  recall  that  equations  (B-l)  through  (B-4) 
came  from  equations  (14)  through  (17),  which  in  steady  state 
per-unit  become: 

(V,)      =  -  —  (A  ) 
d  p.u.      co    q  p.u. 

(V  )  —  (A,) 

q  p.u.      co     d  p.u. 

(V.)      =    JiL  (Xr) 

v  5  p.u.      co     E,   p.u. 

(Vr)      =    —  (A.) 
£  p.u.      co     6  p.u. 

Thus  we  see  that  the  quantities  which  we  must  find  in 
per-unit  are  the  fluxes.   The  flux  equation  is  equation  (21), 
which  is  broken  down  into  individual  flux  equations  below. 

Ad  =  Vs  +  I  M  cos<f)  V  §  M  sin*  h  +  /i  Vf    (B"10) 


(B- 

-6) 

(B- 

-7) 

(B- 

-8) 

(B- 

-9) 
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'Vp.u.    =    (V   /to    )  (B"11) 

r  cr     o 


(X,)  =    ~    [    L„(i,)  I      +    4  M   cos<f>(ix)  I     (N    /N    ) 

d'p.u.         V  E      d  p.u.    o         2  y      6    p.u.    o      a7     a 

-    |  M    sin<|>(ir)  I     (N   /N    )     +     /|  L,  (i-)  (|   I    )  (N    /N,)     ] 

2  Y       £    p.u.    o      a'     a  /    2      1      f    p.u.    2      oaf 

(B-12) 

3  3 

A      =   L^i      +    r  M   sind)    ir    +    t-  M   cosd>    i-  (B-13) 

qEq2  r       &         2  K 

X 

(*     )  =    7TT-7 r  (B-14) 

q'p.u.  (Vo/ooq) 


CO  - 

'Vp.u.    "  v;   ILB(iq>p.u.Io  +   f  M  ■ln*(Vp.u.Io(W 

+  |M  "^'Vp.u-'oW  (B-15) 


Xj    =  i  M  cos*    i  .   +   I-  M   sin<t>    i      +   L   i      +    /!■  L-,i,      (B-16) 
62  a        2  <3  £  /  2      2   f 


(Vp.u.    =      V5    N  (B~17) 

(— )  (— ) 
o        a 


(AP)  =    —  [^Mcoscb  (i  ,)  I     (N    /N    )  +  ^Msinc{)(i    )  I     (N   /N    ) 

v    6    p.u.         V    l2  y      d'p.u.    o       a7     a         2  Y      q    p.u.    o      a'     a 

/T  1  N    2 

+    L     (ij  I     (N    /N    )2+/4L9(i.)^    ln     (yl    )  (;H§-  ^ 

e      6    p.u.    o      a      a         /2    2      f    p.u.    2    o      NfN 

(B-18) 
X      =    -   |  M   sin4>id   +    |  M   cos^i      +    L£i,  (B-19) 
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X 

(VP.u.    =      V        N  <B-20> 

(— )  (— ) 
o        a 


00  -. 

o     r         3 


(Xr)r.    „         =   w~    t    "    T  M   sincf)(i,)  I     (N    /N    ) 

t,   p.u.  V       L         2  Ydp.u.    oa'ct 


+   4  M   sin<j)(i    )  I    (N  /N    ) 

2  Y      q '  p .  u .    o      a'     a 


+   L     (ir)  I     (N    /N    ) 2]  (B-21) 

e£p.u.    oa'a  ' 


Now    identify    the    following   quantities 


/       \  w   L_I 

(x  i)                     o   E   o  ,„    __% 

d'p.u.    =    — (B-22) 

o 


^    co   MI         N 

(xm'P.u.  =  <§^><ir>  <s-"> 

^  o  a 


oo   L    I         N 

<Vp.U.    "     <^V^)<tT>2  (B-24) 

^  o  a 


->        ,_-  oo   L,  I         N 

<xli>p.u.  =  2  /I  ^v^  «n!>  <b-25> 

V  of 


-s        ,-=•   oo   L?I  N2 

(xL2}p.u.    =    2    /    2    ~~ V  (iTN7)  (B_26) 

oaf 


Now    substituting   equation    (B-22)  through    (B-26)    into    the 

foregoing    flux   equations    and    dropping  the   parentheses    and   the 

subscript    "p.u."    since    all   quantities  are   now   in   per-unit 
yields : 
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X,    =   x,i,    +   xi?    cost})   -    x   ir    sinij)   +   xT1i,  (B-27) 

a  a  a  mo  m   £,  LI    f 


X      =   x,i      +   x   i,    sind>   +    x   ir    cosd)  (B-28) 

qdqmom^ 


\,    =   x   i,    cosij)   +   x   i      sin<}>   +   x.i.    +   x_  _i  _  (B-29) 

6  md  Y  mq  r  oo  L2f 


Xr   =-x   i,    sind>    +    x   i      cosd)    +   xri_  (B-30) 

£mdTmqo£ 


Substituting  these  results  into  equations  (B-6)  through 
(B-9)  and  noting  that  for  steady  state  synchronous  operation 
co/oj      =    1.0    p.u.    yields: 


V,    =    -(x,i      +   x   i.    sincj)   +    x   ir    cosd))  (B-31) 

d  d  q  m6  m   t, 


V      =   x,ij    +  x   ix    cosd)    -    x   ir    sind)    +   xT  ,  i,.  (B-32) 

q  d   d  mo  m   t,  LI    f 

Vx    =    x   i,    sind)    -    x   i      cost})    -    xrir  (3-33) 

6mdrmqr6£ 

V^    =    x   i,    cosd)   +    x    i      sin*    +    xrif    +    xT  _i£  (B-34) 

£  md  r  mq  66  L2f 

Similarly,    equation    (B-5)    can  be    put   into   per-unit   as 
follows : 

0    =   —  L, (i    )  I       +    —  L_(ir)  I     (N    /N    )  (B-35) 

co         1      q'p.u.      o         oo         2      £    p.u.       o      a'     a 
o  ^   c  o 

Multiplication  of  both  terms  by  the  proper  constant  and 
assuming  steady  state  operation  yields: 

0  ■  "li1,  +  "l^s  (b-36) 
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Equations  (B-31)  through  (B-34)  and  (B-36)  are  the 
governing  equations  for  the  synchronous  transformer  in  per- 
unit. 

The  following  definitions  are  introduced: 

Vl  -  (vtl>p.u.  <B-37) 

V2  "  <Vt2)p.u.  (B"38) 

Substitution  of  the  definitions  of  equations  (26)  through 
(29)  ,  (B-37)  ,  and  (B-38)  ,  directly  yields  the  per-uni,t  equa- 
tions for  the  synchronous  transformer  operating  between  two 
buses . 


V,sin6      =    x,i      +    x   ir    sin6   +   x    i_    cos*  (B-39) 
ladqm6rmC 

Vncos6      =   x,i,    +    x   i.    coscj)   -   x   i_    sind>  +   xT  .  i  _    ,_    ... 

1           add           mo                        m   £  Llf     (B-40) 


V„sin(6    +5    )    =    -x   i,    sintp    +    x   i      costj)   +    x,ir  (B-41) 

2  a      o  md  Y  mq  6    E, 


V^cos  (5    +6    )    =   x   i,    coscf)   +    x  i      sine})   +   x,i£    +   xT  -i,.    (B-42) 
2  a      o  md  T  mq  o    6  L2f 


0    "    "L^q    +    *L2H  (B"43> 
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APPENDIX  C 
GAUSS-SEIDEL  METHOD 
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GAUSS-SEIDEL  METHOD 


The  Gauss-Seidel  method  is  a  useful  technique  for 
solving  the  following  system  of  simultaneous  linear  equations 


11  1     12  2  In  n    l,n+l 


a^-,xn     i     a  ~  ~  x  _    ■+■    •••    +    a^    x      —    a  «       ... 
21    1  22    2  2n   n  2,n+l 


a    nX,    +   a    n x ~    ■+■    •  •  •         a_    x      —   a  ,  ,      , 

nl    1  n2    2  nn  n  n,n+l  (C-l 

In  order  to  start  the  algorithm  an  initial  solution 
vector  must  be  assumed  and  the  interative  algorithm  is: 

a .   ,  ,     n   a .  . 
x.  =  -i*S+l  -   [   -il  x.  (C-2) 

1       aii         jii  aii     ^ 

In  the  iterative  algorithm,  (C-2) ,  the  most  recently 
available  values  of  x.  are  always  substituted  in  on  the 
right-hand  side. 

The  convergence  criterion  is: 

|x. (k+lst  iteration)  -  x. (kth  iteration) |  <  £         (C-3) 

i  =  1 ,  2 ,    .  •  •  /  n 

For  computer  solution  an  upper  limit  on  the  number  of 
iterations  should  be  specified  since  convergence  may  not 
occur. 
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